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Abstract 

We study the Susceptible- Infected-Recovered model in complex networks, 
considering that not all individuals in the population interact in the same 
way between them. This heterogeneity between contacts is modeled by a 
continuous disorder. In our model the disorder represents the contact time 
or the closeness between individuals. We find that the duration time of an 
epidemic has a crossover with the system size, from a power law regime to a 
logarithmic regime depending on the transmissibility related to the strength 
of the disorder. Using percolation theory, we find that the duration of the 
epidemic scales as the average length of the branches of the infection. Our 
theoretical findings, supported by simulations, explains the crossover between 
the two regimes. 
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1. Introduction 

Complex Networks has became a topic of interest among scientists in 
the last years, due to the fact that many real systems such as protein in- 
teraction, Internet, communication systems [H, 0, S], among others, can be 
properly described by complex networks, making this theoretical framework 
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inherently interdisciplinary. On complex networks, nodes can represent the 
individuals of a population in the case of social networks, computers in com- 
munication systems, etc, and the links represent the interaction between 
them. The research on networks goes from the study of its topology to the 
study of transport processes that use the networks as the underlying sub- 
strate to propagate. In particular, many researchers have focused on the 
propagation of seasonal diseases on social networks due to the appearance of 
new infections such as SARS, Avian Flu and the recently influenza A(H1N1) 
epidemic. 

Several models have been developed to characterize the spreading of this 
kind of diseases, one of the most used models is the Susceptible-Infected- 
Recovered (SIR) model, first introduced by Kermack and McKendrick jij 
in the full mixing approximation and then extended to complex networks 
(5), 0] . In this model the individuals of the population can be in three different 
states S (susceptible) - I (Infected) - R (recovered or removed). A S individual 
becomes infected with a probability by contact with an infected individual. 
Infected individuals recover after t r times steps since they were infected, and 
cannot infect or change their state thereafter. The system reaches the steady 
state when all the infected individuals recover. In this model, the size of the 
disease defined as the number of recovered individuals, depends only on the 
effective probability of transmission between individuals, given by 



where P(/3) is the density distribution of (3, with P((3) = 8{j5 — for a 
constant probability of infection. 

It has been shown that the steady state of the SIR model on statics 
complex networks can be mapped into a link percolation problem js), 0] where 
the order parameter of the SIR is the fraction of recovered individuals, and 
the control parameter is the transmissibility T which plays the role of the 
fraction of links p in percolation. Thus there exist a critical threshold T c 
(or p c ) in the SIR (or percolation) model above which a nonzero fraction of 
individuals (size of biggest cluster) are infected (~ N). It was shown that 
in finite systems, the critical threshold T c = p c of link percolation depends 
strongly on the network size N with T C (N) — T C (N — > 00) ~ iV -1 / 3 in mean 
field approximation, i.e. finite size effects are strong [8j. Then in finite 
systems, for T > T C (N) the disease becomes an epidemic, while for T < 
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N) the disease reaches a small fraction of the population (outbreaks) 



a 10, 11 



Usually, the study of a disease transmission or any type of transport pro- 



cess like information flows 12 



13| through a network is 
However this assump- 



rumor spreading 
made assuming that all the contacts are equivalent 
tion does not give a very realistic description of real networks, such as social 
networks 14( where not all individuals in a society have the same interaction 
between them. A way to improve the description of real networks, is to con- 
sider the heterogeneity of the social contacts. This can be done by considering 
weighted (disordered) networks, where the weights in the contacts could rep- 



resent the closeness or the contact time between the individuals [15|, [16|, 111 
The contact time is a parameter that can be controlled by health policies as 
a strategy to mitigate the duration of the spreading of the disease. Using 
different strategies like broadcasting, brochures, etc, the public health agen- 
cies can induce people to change their contact time or the closeness of the 
contact, for example, encouraging people to reduce their contact time. This 
strategy that is a social distancing, was used already by some governments 
in the recent wave of influenza A(H1N1) epidemic in 2009 



18 . It is known 



that the disorder can dramatically alter some topological properties of net 



works such as the average length of the optimal paths [15|, [19|, |20|, |2l| . In 



the optimal path problem defined as the path between any two nodes that 



minimizes the total weight along the path [15l . Il9l . |20| , it was shown that the 



average length of the optimal path l opt scales as N Uopt in the strong disorder 
limit (SD), where z/ opt = 1/3 for homogeneous networks, and as l op t ~ IniV in 
the weak disorder regime (WD), where the SD limit is related to percolation 



at criticality |15|, [19|, |22j. However, the exact mapping between the order pa- 
rameter of both second order phase transitions, percolation and SIR, is not 
affected by the disorder when the disorder is not correlated (§]. Nevertheless, 
the disorder could affect the duration of a disease. 

In this paper we introduce disorder in the links and show how a broad 
disorder affects the duration of an epidemic in the SIR model. 

Using theoretical arguments, supported by intensive simulations, we find 
that, the average time of the duration of the epidemic tf goes as the average 
length of the branches lb of the infection. Thus, relating Ij, with the optimal 
path problem we find that tf has a crossover from WD to SD regimes. 

The paper is organized as follows: In section [2] we present our model, 
in section [3] we show theoretical and by simulation, the crossover for the 
duration time of the disease from WD to SD regimes. Finally, in section H] 
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we present our conclusions. 



2. SIR model with disorder 

In our model, we assign to each link i between any two nodes a random 
number drawn from the distribution 

p(Pi) = 4r> ( 2 ) 

aPi 

where a is the parameter that controls the broadness of the distribution of 
link weights i.e. the strength of the disorder and is defined in the interval 
[e~ a , 1]. With this distribution we randomly assign a weight to each link of 
the network, of the form, 

Pi = e~ ar \ (3) 

where r« is a random number taken from an uniform distribution ri e [0, 1]. 
Thus fa is the probability of infection between any pair of nodes. This type 



of weight has been widely used [19|, |22|, |2l|, |23j because it is a well known 



example of many distributions that produce WD and SD crossover. The 
strong disorder limit can be thought in a disordered medium as a potential 
barrier q so that T{ is the time to cross this barrier in a thermal activation 
process, then n = e 6i ' kT , where k is the Boltzmann constant and T is the 
absolute temperature. Therefore, we can see as the inverse of the time 
needed to cross this barrier Tj. 

With this weight distribution the transmissibility T (See Eq. CP is given 

by, 

In our initial stage all the individuals are in the susceptible state. We choose 
a node at random from the biggest connected cluster, or giant component 
(GC) and infect it (patient zero), then the process follows the rules of the 
SIR model but in a weighted network. After the system reaches the steady 
state, we compute the duration of the infection t* defined as the time where 
the last infected individual recovers and the length of the branches of the 
infection as a function of t*. As the substrate for the disease spreading 
we use only an Erdos Renyi network (ER) j24j|, characterized by a Poisson 
degree distribution P(k) = (k) k k\ where k is the connectivity and (k) is 
the average degree. However, the results obtained are similar for all networks 
with a finite value of T c in the mean field approach. 
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Figure 1: P(t*) as a function of t* for an ER network with (k) = 2, N = 2 14 , 
t r = 19, a = 4.36 (full line) and a = 3.24 (dashed line). The arrows indicate 
the separation between outbreaks (left) and epidemics (right) regime. In the 
inset we show in symbols the fitting with Eq.Q of the data for a = 3.24 (see 
main plot) from where we obtain t/ — 363.8. All simulations were done over 
10 5 realizations. 



3. Crossover from WD to SD for the duration times of the epi- 
demics 

We only consider those propagation that lead to epidemic states, and dis- 
regarded the outbreaks (See Fig. [Q. Fig. [1] shows P(t*) as a function of t*, 
the arrows shows the separation between outbreaks and epidemic regimes. 
The criteria used to distinguish between outbreaks and epidemics is the obser- 
vation of the behavior of the cluster size distribution that decays as a power 
law close to T c for the outbreaks and has a maximum for the epidemics. The 
same criteria is used in percolation to distinguish finite cluster sizes from 
giant component with size ~ N. We find that P(t*) in the epidemic regime 
can be well represented by a Log-Normal distribution, 



P(t* 



1 



ln(t*/t f ) 2 



t*aV2TT 



(5) 



where tf and a are the average and the standard deviation of the distribution 
of t*. Log-normal distributions have been observed in several phenomena, 
such as the size of crushed ore 25(, fragmentation of glass 26(, income dis- 



tribution 27 



mastication 



29 



events in medical histories 28], food fragmentation by human 
etc. By fitting our data with a log- normal distribution we 



obtain tf (see the inset of Fig. [T]). 
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Figure 2: Log- log plot of the integral of 4 as function of t*, for an ER network 
with (k) = 2, N = 2 14 and t r = 19. For the undisordered problem T = T c 
(□), and for the disordered problem with T = T c and a = 6.55 (0)> an d 
T > T c and a = 3.97 (A). The dashed line is used as a guide to show the 
slope 2. All simulations were done over 10 5 realizations. (Color online) 



From Fig. [T]we can see that tf increases with a for a fixed value of t r , this 
behavior can be understood if we take into account that when a increases, 
the barrier that the disease needs to overcome in order to propagate is bigger. 
Therefore, even though the transmissibility decreases as a increases for fixed 
t r , the epidemic will last longer due to the disorder, allowing the health 
services to make earlier interventions. 

We also study the average length of the branches of the infection The 
length of a branch is defined as the number of links between the patient 
zero and the last patient in that branch. In Fig. [2] we plot the integral 
of lb as a function of t* in log-log scale, from which it is easy to see that, 
independent of the disorder, the epidemic region has a slope 2, thus i& ~ t*. 
If we compare lb between the disordered and the undisordered substrate, we 
can see that in the disordered one, the length of the branches are larger than 
in the undisordered case, even for the same values of T and t r . This is due 
to the fact that the link of smallest crossing probability e~ a that the disease 
has to traverse is smaller than the required value of /3 for the same values of 
T and t r in the undisordered case (see Eq.flTJ). 

For the SIR problem, at the threshold T c in the mean field approximation, 
the fraction of recovered individuals R grows with time t* as R ~ t* 2 [30j up 
to tf, then at tf, R ~ t 2 f, and for t* > tf, R ~ constant due to finite size 
effects. Percolation in complex networks predicts that at criticality T w T c , 
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Figure 3: For an ER network with (k) = 2, t r = 19 and T = T C (N = 2 12 ). 
(a) Linear-linear plot of tf/N ' 36 as a function of N, the dashed line is used 
as a guide to show that the exponent 0.36 is correct, (b) Linear-log plot of 
t f as a function of ln(iV), the dashed line correspond to a logarithmic fit. 
Notice that in the WD regime T C (N = 2 12 ) > T C (N > 2 12 ). All simulations 
were done over 10 5 realizations. 



the size R of the epidemics scales with the chemical length I as [31 



with di = 2 in the mean field approximation, then 

t f ~i, 

where t ~ A^ 1/3 . For the WD regime where T > T c 



R 



then at t 



f 



e<* for t*<t f ; 
const, for t* > tf , 

R ~ exp(tf) . 



(6) 



(7) 



(8) 



(9) 



Using the fact that in WD , R ~ iV and i ~ In N, we obtain that also in this 
regime 



(10) 
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Thus from the theoretical arguments presented above, 




N i/3 if t~T c 



(11) 



ln(JV) if T>T C 



In order to corroborate Eq. ffTTl) . we compute tf as a function of N for 
fixed values of T = T C (N = 2 12 ), a and t r , and found that for T ~ T C (N), tf 
behaves as a power law with exponent 0.36 compatible with the one found for 
SD in the optimal path problem. The slight difference between the exponent 
obtained and the theoretical one (1/3) is due to finite size effects |8J. Notice 
that the power law regime of SD holds only until iV ~ 15000. In Fig. [3] (a), 
we plot the ratio tf/N - 36 as a function of N. We can see that for values 
of N for which T ~ T C (N), this ratio goes to a constant in agreement with 
our prediction. However for T > T C (N) we can see from Fig. [3] (b) that tf 
behaves logarithmically with N in agreement with the WD regime T > T c of 
Eq. ffTTT) . Thus, we can see clearly that there is a crossover from WD to SD 
in full agreement with our theoretical prediction. 

4. Conclusions 

In this paper we study the SIR model on a broad disordered network, 
where the disorder represents the duration of the interactions between the 
individuals of a population, and the broadness of the disorder is a control 
parameter. 

We found that as the broadness of the disorder increases, the spreading of 
the disease is delayed. Thus, recommending people to decrease the duration 
of the contact as a social distancing could be a good strategy to delay the 
spreading of an epidemic. Moreover, this social distancing strategy, is more 
suitable than a quarantine, where contacts are broken, due to the fact that 
it is less expensive from an economic point of view. 

Using percolation arguments, we found that the duration time of the epi- 
demic goes as the average length of the branches of the infection. Our theo- 
retical results are strongly supported by simulations. Thus, in the same way 
as in the optimal path problem, the duration of the disease has a crossover 
with the system size, from a power law regime (SD) to a logarithmic regime 



(WD). 
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